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Abstract 

We give a characterisation of the spectral properties of linear differential operators with 
constant coefficients, acting on functions defined on a bounded interval, and determined by 
general linear boundary conditions. The boundary conditions may be such that the resulting 
operator is not selfadjoint. 

We associate the spectral properties of such an operator S with the properties of the 
solution of a corresponding boundary value problem for the partial differential equation dtqi 
iSq = 0. Namely, we are able to establish an explicit correspondence between the properties 
of the family of eigenfunctions of the operator, and in particular whether this family is a basis, 
and the existence and properties of the unique solution of the associated boundary value 
problem. When such a unique solution exists, we consider its representation as a complex 
contour integral that is obtained using a transform method recently proposed by Fokas and 
one of the authors. The analyticity properties of the integrand in this representation are 
crucial for studying the spectral theory of the associated operator. 

MSC: 47A70, 47E05, 35G16, 45P10, 35C10 

1 Introduction 

In this paper, we study the following two objects: 

(1) A linear constant-coefficient differential operator S defined on a domain of the form 
T)(S) — {u £ L 2 [0, 1] : u sufficiently smooth and satisfying n prescribed boundary 
conditions}. 

(2) An initial boundary value problem (IBVP) for the linear evolution partial differential 
equation qt(x,t) ± iSq(x,t) — 0, x G (0,1) t 6 (0, T), with S as in (1), an initial 
condition q(x,0) = qo(x) and n given boundary conditions. 

The boundary conditions, assumed to be linear, can be prescribed at either end of the interval 
[0, 1], or can couple the two ends. 

It is to be expected that the objects (1) and (2) are closely related. For each of these 
objects, it is natural to formulate a basic question, whose answer depends on the specific 
boundary conditions. Namely, given a set of n boundary conditions, 

(Ql) does the resulting operator S admit a basis of eigenfunctions, in any appropriate sense? 

(Q2) does the resulting initial-boundary value problem admit a unique solution representable 
by a discrete series expansion in the eigenfunctions of S? 
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Although it should be clear that these are the same question posed in different contexts, 
very little is explicitly known beyond the classical cases when the spatial operator has a 
known basis of eigenfunctions. This basis can be used after separation of variables to express 
the solution of the boundary value problem. 

In this paper we give an explicit connection between the two problems in general; we give 
a link between the solutions of (1) and (2), and we show precisely how the answer to (Ql) 
and (Q2) are related. In particular, the rigorous answer to one question can be given through 
answering the other. Our results are true for general n, however they are new and interesting 
in particular for n odd. 

Since in general S will not be self-adjoint, we expect that any spectral decomposition 
involves not only S but also the adjoint S* . In terms of the PDE problem, we will see that 
this is reflected in the need to consider both the initial time and the final time problems (the 
evolution with reversed time direction). 



The operator problem 

We consider the linear ordinary differential operator S, given by 

d 



Su = 



ax , 



uev{S), 



defined on the domain T>(S) C L [0, 1] given by 

V(S) = \u G AC n [0, 1] : A(u (B_1) (0), ^""''(l), . • ■ , u(0),u(l)) = 0}, 

where 



(1.1) 



(1.2) 



AC n = {/ G C" 1 " 1 : / (n_1) absolutely continuous and / (n) G L 2 [0, 1]}. (1.3) 

By T>(S), we denote the L 2 closure of T>(S). Here the order n ^ 2 is an integer and the 
boundary coefficient matrix A G R" x2n , encoding the given boundary conditions, is of rank n 
and given, in reduced row-echelon form, by 



A : 



/ain-1 Pi n — 1 Q'ln-2 Pin — 2 
OL2n-l $2 n — 1 0^2 n — 2 P2 n — 2 



Pnn—1 Ct nn -2 Pnn—2 



aio Pio\ 

OL20 P20 
OtnO PnO/ 



(1.4) 



coefficients. 



The numbers ctj r , Pjr are called the boun 

This operator has been studied at least since Birkhoff (1908b). Depending on the particular 
entries of the matrix A, the operator may or may not be selfadjoint. The theory of the 
sel fadjoint case was ful ly understood by the time lDunford and Schwartzl (| 19631 ) presented it. 

I Locked (|2000l . I2008T ) used the theory of Fredholm operators to study the non-selfadjoint 
case. He defined the characteristic determinant 

A(p) = e l ^=' u ' kp detM(p), (1.5) 

where oj = exp(27ri/n) and the entries of the matrix M(p) are given by 

71 n 

M kj{p) = ^ajr(itJ fc_1 p) r e" 1 " " + ^i3 jr {iu k ~ 1 p) r ■ 

r — 1 r — 1 

It is known that, provided A 7^ 0, if A(er) = then a n is an eigenvalue of S. Further, 
the algebraic multiplicity of a n as an eigenvalue of S is equal to the order of a as a zero of 
A. Locker showed that, for Birkhoff-regular operators, the generalised eigenfunctions form a 
complete system. However, he gives no general statement about the cases that do not satisfy 
these regularity conditions. 



2 



The PDE problem 



In a separate development, a novel transform method for analysing IBVPs was developed 
by Fokas (see Fokas, 20qI, for an overvi ew). The method wa s applied to IBVPs posed 
for evolution equatio n s on the half-line bv lFokas and Sunel (| 1999h and on the finite i nterva l 
by iFokas and Pellonil (|200lh with simple, uncoupled boundary conditions. In ISmithl |2012), 
Fokas' method was applied to IBVPs whose spatial part is given by the operator S, namely 
those of the form 

d t q(x,t) + a(-id x ) n q{x,t) = 0, x G (0, 1), t > 0, a = ±i, (1.6) 

with prescribed boundary conditions and an initial condition qo(x) = q(x,0), assumed smooth 
to avoid technical complications. Usually the initial condition is assumed to be in C°°. How- 
ever, the same results hold assuming that go G AC n . Indeed, in this case, the uniform 
convergence of the integral representation (see (|1.7p below), the poynomial decay rate of the 
integrand and the explicit exponential x dependence imply that the solution q belongs to the 
same class. In what follows we assume go G AC n . 

This method yields an integral representation of the solution of the initial-boundary value 
problem in the form 

q{x,t) = ^-[ e^- ap " i 4^dp+— / e 1 *— ^--'^-p^-dp 

2-7T J r+ ApDE 27T J r - ApDE 



+ i *y e i°-k*-»"?« Res ^S£l dp + i -y e io-*(*-i)-a*?t Res p^-dp, (i.7) 

P="h ApDE Z— ' P=°k ApDE 



where the quantities go, C , Apde , Ofc and Y are defined below in Definitions 12. II and 12.41 In 



many c ases, inc luding all problems with n even, the integrals in equation (|1.7[) both evaluate 
to zer o (ISmith L 2012)'). W e study these cases here. 



In lPellonil if2004U2005l ) and then in greater generality in lSmitr] (|201ll ). this method is used 
to characterise boundary conditions that determine well-posed problems, and problems whose 
solutions admit representation by series. To achieve this characterisation, the central objects 
of interest are the PDE characteristic matrix A (see Definition 12 . 1 1 below) and its determinant 

ApDE • 

Note that in this work, by 'well-posed', we mean existence and uniqueness of a solution 
and make no claim to continuit y with respect to data. By 'ill-posed' w e mean that ex istence or 
uniqueness fails. The results of lFokas and Sund jl999h : IPellonl l|2004h : ISmitnl (|2012h establish 
that a problem is well-posed if and only if it admits a solution via the method of Fokas. 

The present work details results connecting the spectral theory of S with the behaviour 
of the associated IBVPs for the PDE H1.6[) . as well as the one obtained from the same set of 
boundary conditions but for the PDE 

d t q(x, t) - a{-id x ) n q{x, t) = 0, x G (0, 1), t > 0. (1.8) 

We refer to the latter in the sequel as the final time boundary value problem. 



Summary of the main results 

For an operator S of the type given by (|1.1[) . and the associated initial- and final-boundary 
value problems, we prove the following: 

• If the eigenfunctions of S and S* form a biorthogonal basis of T>(S) and the IBVP is 
well posed, then its solution is representable as a series. 

This is the content of Proposition ^. 71 It follows from this result that if a series represen- 
tation does not exist, then the eigenfunctions of S and S* cannot form a basis of T>(S). 
What is interesting is that we can use the PDE approach to obtain results on S in cases 
that are not covered by usual operator theoretic techniques. In section 0] we provide an 
example when (Ql) cannot be answered by the usual tests involving projector norms, 
but may be settled through this result and a negative answer to (Q2). 
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• If the initial- and final-boundary value problems are well posed, then the eigenfunctions 
of S and S* form a complete biorthogonal system in T>(S). 

This is the content of Theorem 12.61 The conclusion does not imply that the eigenfunc- 
tions necessarily form a basis. However the integral representation (|l,7|l can always be 
deformed to derive a series representation for the solution of the IBVP in terms of the 
eigenfunctions. 

• The departure of the family of eigenfunctions of S and S* from being a biorthogonal basis 
can be estimated in terms of the integrand in the representation of the solution of the 
associated IBVP. 

This is the co ntent of Theo rem 12.121 This departure is quantified in the notion of 
'wildness' (see iDaviesl 120071) . Indeed, if the eigenfunction of S and S* form a wild 
system in L 2 [0, 1], then we provide an estimate of the wildness of the system in terms of 
the quantities used to determine whether the initial- and final-boundary value problems 
are well posed. 

Outline of paper 

In section [2l we review the necessary definitions and notation. Following this, we precisely 
state and prove the results described above. 

Each of sections [3] and U is devoted to the analysis of an example which illustrates the above 
general results. We compare and contrast the results obtained through the new theorems with 
those yielded by Davies' wildness method. 

2 Complete and basic systems of eigenfunctions 
2.1 Notation, definitions and preliminary results 

In this paper, we make extensive use of the notation developed in ISmithl (|2012l ). We refer to 
that paper for details, but we list here some of the notation used throughout the rest of this 
work. 

The initial-boundary value problem U(n,A,a,q ): Find q G AC n ([0, 1] x [0,T]) which 
satisfies the linear, evolution, constant-coefficient partial differential equation 

d t q(x, t) + a(-id x ) n q(x, t) = (2.1) 

subject to the initial condition 

q(x,0)=q (x) (2.2) 

and the boundary conditions 

a (ar Vo, t), ar Vi, t), ar 2 g(o, t), ar Vi, *)>•--, a(o, t), *)) T = Kt), (2.3) 

where the quadruple (n,A,a,q ) G N x R nx2rl x C x AC n [0, 1] is such that 

(111) the order n > 2, 

(112) the boundary coefficient matrix A is in reduced row-echelon form, 

(113) the direction coefficient has the specific value a = ±i, 

(114) the initial datum qo is compatible with the boundary conditions in the sense 

A (cfr- 1} (0), q^ x) (1), 9 «- 2) (0), ^- 2) (1), . . . , go(0), g (l)) T = 0. (2.4) 

Given a problem II = II(n, A, a, qo), we define the corresponding final time time problem 
IT = IL(n,A,-a,qo). 

We assume that the boundary conditions are homogeneous to aid the comparison with 
S, the differential operator representing the spatial part of the PDE problem II. There is 
no loss of generality in this assumption. Without this restriction, II is no more difficult to 
solve ; the solution simply contains an additional term represented as an integral along the real 
line l|Smithl . l2012h . 
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Definition 2.1. Let a\p ft^i be the boundary coefficients of the operator S* , adjoint to S. 
We define 



4W = E(- #1 f)^i. ( 2 - 5 ) 

r = 
n-1 

^■(^Ef-^Vffi, (2.6) 

then A kj (p) = Aifr) + A^ j (p)e~ iuk ~ lp (2.7) 

is called the PDE characteristic matrix. The determinant Apde of A is called the PDE 
characteristic determinant. 

Remark 2.2. The PDE characteristic matrix is a realisation of Birkhoff 's characteristic ma- 
trix for S* and also represents the Dirichlet-to-Neumann map for the problem EL Indeed, 
it is through this matrix that the unknown (Ne umann) bound ary values are obtained from 
the (Dirichlet) boundary data of the problem. ISmithI ||2012T ) uses a different but equiva- 
lent definition of A which generalises the construction via determinants and Cramer's rule 
originally found in 



m 



Fokas and Sund i|l999l ). The validity of the new definition is es tablished 



Fokas and Smith! (I2013f 7 and the equivalence is explicitly proven in lSmithl l|2013bf ). 



Remark 2.3. In Definition 12.11 we construct A via the boundary conditions of S*. It is 
possible to make an alternative but equivalent definition of A via an explicit construction 
from the boundary conditions of S itself. For the examples considered in sections [3H3] this 
i s a s imple matter. Indeed, provided the boundary conditions of S are non-Robin, ISmithI 
(|201ll . Lem ma 2.14) pro vides a simple construction. This can be done for general boundary 
conditions (Smith, 2012) and can easily be coded to be done automatically. 

Definition 2.4. Let (o~k)kew be a sequence containing each nonzero zero of Apde precisely 
once. We define the index sets K + = {k G N : a k G C+}, K~ = {k G N : o k G C~}. Let 
3e be the infimal separation of the zeros a k . Then the contours T are the positively-oriented 
boundaries of 

{pGC* : Re(ap n ) > 0} \ \J B(a k ,e). (2.8) 

The minor X r 3 (p) is the (n— 1) x (n— 1) submatrtx of A whose (1,1) entry is the (r+l,j+l) 
entry. This is used to construct the spectral functions 

n n 

r=l 3=1 
n n 

C(p,qo) = Y,Y / detXrj (P') A Up)Q°(" r ~ 1 p)> ( 2 - 10 ) 

r — 1 j = l 



where 



l 

ipx 



qo(p) = / e lpx q (x)dx. 



(i 



Definition 2.5. We say the IBVP is well-conditioned if it satisfies: 
C ± (p) is entire and the ratio 

f ± (p) as p — > oo from within a sector exterior , , 



. (p) to r , away from the zeros of Apde • 



Otherwise, we say that the problem is ill-conditioned. 

Well-conditioning of an IBVP is not a classical definition and is unrelated to the concept of 
conditioning that appears in numerical analysis. Conditioning, in the sense of Definition 12.51 
is necessary for well-posedness but is also central to the validity of a series representation. 
Indeed, switching the direction coefficient a h-> —a in the PDE (|1.6p switches which sectors are 
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enclosed by the contours T ± thus, by Jordan's Lemma, well-conditioning of the proble m with 
the o pposite direction coefficient is equivalent to the two integrals in (|1.7[) vanishing l|Smithl 
I2012D . 

The reader will recall that a system ((d„)„ 6 i in a Banach space is said to be complete if its 
linear span is dense in the space and such a system is a basis if for each / in the space there 
exists a unique sequence of scalars (a„)„ 6 N such that 



/= lim (£\ 

\n=l 



2.2 Well-posed PDE systems and bases of eigenfunctions 

It is well known (see ICoddington and Levinsonl 1 19551 . Section 12.5) that if the zeros of the 
characteristic determinant A of S are all simple then the eigenfunctions of S form a complete 
system in T>(S). This theorem is proven using an analysis of the Green's functions of both the 
operator S and its adjoint S* . We prove the following result without directly analysing the 
adjoint operator. 

Theorem 2.6. Let S be such that the zeros of Apde are all simple. Let II = II(n, a, A, qo, 0) 
be an IB VP associated with A and II' be the corresponding problem with the opposite direction 
coefficient, II(n, — a, A, qo, 0). //II is well-posed and II' is well-conditioned in the sense of 
Definition \2.5\ then the eigenfunctions of S form a complete system in T>(S). 

Rather than analysing both the original operator S and the adjoint operator S* , one 
needs information on both the initial- and final-boundary value problems associated with the 
operator S. 

A stronger, but essentially straightforward, result in the reverse direction is: 

Proposition 2.7. // the eigenfunctions of S form a basis in T>(S) and, for some a, the 
associated IBVP II is well-posed, then II' is well-conditioned. 

Further, if (<f>k)k£N are the eigenfunctions of S, with corresponding eigenvalues (a^)kefi 
then there exists a sequence (ipk)k£?i biorthogonal to (4>k)ksm such that the Fourier expansion 

^•MzXgo^e-^ (2.12) 

feeN 

converges to the solution of II. 

Indeed, in the notation of Proposition 12.71 each ij)k is an e igenfunction of the adjoint 
operator S* with corresponding eigenvalue — cr^ ( Birkhofii. 1 1908al ) . 

The a bove results a re essentially the translation into operator theory language of results 
proved in ISmithl 4 201 il l. Here we extend the parallelism between PDE and operator the- 
ory in important ways. Namely, under some further assumptions, we construct explicitly 
the eigenfunctions of the differential operator directly from the PDE characteristic matrix. 
The construction does not require knowledge of the integral representation even implicitly, as 
neither II nor II' need be well-posed. 

In the sequel, we assume that the boundary conditions are non-Robin and that a technical 
symmetry condition always holds, see Conditions I A. II and l A. 21 in the appendix. We also define 

n 

Q(p;qo) = £detjr i (p)ft)(w p - 1 p), (213) 

r=l 

so that 

n 

C ± (p;go) = ^- 4 ^(/ 9 Ki(p;?o). (2.14) 

3=1 

In the next proposition, we characterise the eigenfunctions of S in terms of the PDE 
characteristic matrix and the spectral functions. 
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Proposition 2.8. For each k £ N and for each j € {1, 2, . . . , n} ; the function 

n 

= J2 e ~^ r ^ k(1 ~ x)detxrj ( ak ) ( 2 - 15 ) 

r = l 

is an eigenfunction of S with eigenvalue er£ . Further, 

= 7ri<lo,ipi}, j = l,...,n, fceN (2.16) 

G(^,go) = Q(g ,^>, (2.17) 



V£(l-x) = C^i(x), (2.18) 

where is the corresponding eigenfunction from the adjoint operator and Cj is a nonzero 
real scalar guantity depending only upon j. 

Remark 2.9. The proposition above requires that the boundary conditions be non-Robin 
and obey the symmetry condition. These requirements may not be sharp but we have been 
unable to find an example failing either condition for which the result holds. 

By Proposition 12.81 the spectral functions of the original and adjoint problems, which we 
denote by Q, obey the identity 

The function go (x) denotes the initial datum of the IBVP. Hence it can be chosen arbitrarily 
in T>(S). The particular choice qo(x) = 4>i( x ) * s admissible since ip^ipo) is C°° by definition. 
With this choice, equation (|2.19p yields 



CjK,^) _ \m\ 2 _ \CAUU\HI\ 



= |C J |||Q fc ||, (2.20) 



c;(^,^ J J {HAD 

where is the projection operator 

QkU) = UAk)i> k (2.21) 

considered bv lDaviesl (|2007l ). Note that the latter equality follows from equation (|2.18p . 
By a simple change of variables we find 



C*(p, *,(■)) = -CMPM^ - •))• (2-22) 

We therefore deduce the following important result, which gives a way to control the norms of 
the projection operators Qk explicitly in terms of the spectral functions associated with the 
corresponding initial and boundary value problem. 

Proposition 2.10. Let S be the operator associated with II. Then the eigenfunctions <j>k and 
tpi of S and of its adjoint satisfy 

II II 2 = CjQ(<rk,ipi) (2 23) 

Remark 2.11. This result implies that we can estimate ||Qfc|| using only the spectral functions 
of the initial- and final-BVPs, whose construction is algorithmic. 

Conversely, this proposition has an important consequence, namely an estimate on the 
unboundedness of the spectral fu nctions i n term s of the "wildness" of the family of biorthogonal 
eigenfunctions of S. (Following Davie^l (|2000l ). we say that a biorthogonal system is wild if 
the corresponding projection operators are not uniformly bounded in norm.) We illustrate the 
result of this theorem in the two examples we consider in sections 3 and 4. 

Theorem 2.12. Let go be any admissible initial condition for the boundary value problem, 
and let {pk)k&t be any sequence such that 
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• pk — > oo as k — > oo. 

• \Pk\ < \pk+i\ 

• (pfc) is bounded away from the set of zeros o/Apde, uniformly in k: 

35 > : V k,j e N, \\p k \ - \aj\\ > 5 



Then 



\\Qk\\ = o( sup 
i (pfc) 



A r 



A r 



as k 



2.3 Sketch of proofs 

Proof of Theorem\2Jk As II is well-posed and II' is well-conditioned, bv lSmithl l|2012l . l2013ah 
the solution q of the problem II can be expressed using a series as 



q(x, t) — i \ Res 



k£K+ 



P = cr k Af 



:(P) 



C + (p) + * Y\ Res 

^ * /l— ITp 



ip(x-l)-ap™i 



-r(p)- 



As each a*, is a simple zero of Apde , the series is separable into a;-dependent and f-dependent 
parts 



Tk(t) = • 



: Res p = c 



C+(p) 

k ApDE (p) 

C~(P) 



so that 



2 e ives p=(Tfc ApDE(p) 



iike K+, 
lik e K~, 



(2.24) 
(2.25) 

(2.26) 



Further, ISmithl l|20ll Lemma 6.1) guarantees the existence o f a no nzero complex constant C 
such that <T fe = Ck + O(l) as k -s- oo, which, by ISedletskiH |2005l . Theorems 3.3.3 & 4.1.1), 
guarantees that (rk)kef> is a minimal system in L 2 [0,T]. 
As q is the solution of II, q satisfies 



A 



/^- 1 ?(0,t)\ 
fl& -1 «(l,*) 



?(0,t) 
V 9(1, *) / 



= 0, 



Vt G [0,T]. 



The minimality of the f-dependent system means that this implies each satisfies the bound- 
ary conditions of S, so £fc G T3(S). 
As g satisfies the PDE, 



= a^[-^7 + S](a)(x)r fc (t) 



by minimality of (Tfc)kgN, each is an eigenfunction of S with eigenvalue cr^ . 
Evaluating equation (|2.26l) at t = yields an expansion of go in the system (£fc)fcgi* 



□ 



Remark 2.13. We have to require the zeros of Apde are all simple. It would be desirable to 
be able to say that the zeros of A and Apde are all the same and of the same order. It h as been 
show n that this holds under certain symmetry restrictions on the boundary conditions l|Smithl 
l201ll ) and has been established in particular for all possible 3 rd order boundary conditions. 
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Proof of Proposition 



As (4>k)k£f>i is a basis, the Fourier expansion 



converges. By I Smith! |2013aT ). well-posedness of IT guarantees that the are arranged in 
such a way that the exponential functions e £ " Tfct are bounded uniformly in k, hence that the 
series (|2.12|) converges for all t € [0, T]. The eigenfunctions all satisfy the boundary conditions 
of the operator so the Fourier series satisfies the boundary conditions of the initial-value 
problem. The Fourier series also satisfies the partial differential equation. So we have a series 
representation of the solution and IT must be well-conditioned. □ 

Proof of Provosition \2.8[ Let Bi be the I th boundary condition of S. As the boundary condi- 
tions are non- Robin, they each have an order mi. Hence 



n 

j) = ^(iafca/- 1 )" 1 ' [a imi e^"~ lCT * +# m ,] detX r] (a k ). (2.27) 



The bracketed expression is an entry from row r of the characteristic matrix of S. Provided 
the boundary conditions also satisfy the symmetry condition, an algebraic manipulation yields 
that each c olumn of the characteristic matrix of S is a scalar multiple of a column of A (see 
the proof of lSmithllioTl] . Theorem 4.15). So either Bi((j>{) is the determinant of a matrix with 
a repeated column or Z?;(<^.) = Apde (o"*0- ^ n either case, Bi(<p 3 k ) = 0, so <j>k £ "D(S). Finally, 

Let the map r n> f be given by the permutation (l,n,n — 1,...,3, 2), whose sign is 
(— l)r n / 2 l _1 . Because the boundary conditions obey Conditions I A . HtA .21 



Cj det X J (p) = e" 
where the real constant 



~'detX rj *(p), VpGC 



(2.28) 



fn/21-l 



l±3 



and Pi is the coupling constant appearing in the I column of A (1 if there is no coupling 
constant in that column). 

Indeed, as S is a closed operator, densely-defined on L 2 [0, 1], the eigenval ues of S* are 
the points a£, and effc are the zeros of the adjoint PDE characteristic matrix i|Smithl. [ioTll . 
Theorem 4.15). Note also that the construc tion of the adjoint boundary con ditions from the 
boundary conditions of the original problem (jCoddington and Levinsonlll955| . Theorem 3.2.4) 
ensures that the adjoint boundary conditions also satisfy Conditions I A . lf[A .21 

As the boundary conditions are non-Robin, the only columns that may appear in A are 



,(n-l)i 



V 



\ 



( 



(e- ip + /3 m ) 

Xn—Vjl r -iu n ~ 



\ 



(2.29) 



where / may vary over {0, 1, . . . , n — 1}. To each of these corresponds a unique column in A* 
with the same values of I as each column in A: respectively, 



( 



\ 



(e- l " + l//3 m ) \ 
w l (e- i *"> + l/p m ) 



,(«-!)! 



(e- 



(2.30) 



Hence, to construct A(p) from A*(p) we apply the following operations: 



1. For all r, multiply the r* row by e u 
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2. For all m, multiply the m th column by /3 m . 

3. Apply the permutation ri->rto the row index. 

4. Take the complex conjugate of each entry. 



This justifies equation (|2.28[) . 

By equation (|2.15[) . the eigenfunctions of the adjoint operator are 

n 

ipi(x)=j2 e ~ lhjrlWk<1 ~ x)detxr] *(° k y ( 2 - 31 ) 

r=l 

By the definition of r t-¥ r, co 1 ~ r — u r ~ 1 . Hence 

n 

ii(x)=^2e i(u ' 1 ~ rx -^~ 1)w '°detX*' j *(a k ). (2.32) 

r=l 



Hence, by equation (12. 2811 . 

n 

= C^e^^^^detX^K) 



= C J^2 e ~^ T 1<TkX detX r 3(a k ). 
Hence, by the definition of it follows that Cj(°"fc) = {lo,tpi)/Cj. □ 



3 Third order coupled and uncoupled examples 

In this section we outline the analysis of a particular class of boundary value problems, de- 
pending on a real parameter ft, for the third order PDE q t = q xxx . Namely we consider the 
following problem: 

qt = q xxx , 35 e [0,1], te[0,T], (3.1) 

q(x,0) = qo(x), x G [0, 1] 

q(0,t) = q(l,t) = 0, q x (0,t)+Pq x (l,t) = 0, t € \0,T\, G R 

where go G D(S) is a known function. 

In the limit as the constant /3 — >■ 0, the second boundary condition at x — is q x (0, t) = 0. 
The spectral properties of this limiting case are very different from the case /3 ^ 0, when the 
coupling between the first order derivatives is lost. Hence we refer to the boundary conditions 
corresponding to the value /3 = as uncoupled. 

In this section we analyse the behaviour of the associated differential operator in the two 
cases. To avoid technicalities, and to concentrate on the ft — limit, we assume in what 
follows that /3 G (-1, 1). 



The associated differential operator 

Let S 13 be the differential operator corresponding to the boundary value problem ()3.1|l , hence 
specified by n — 3 and by the boundary coefficient matrix 

/0 1 p 0\ 
A p = 1 o , pe (-1,1). (3.2) 



Setting 



v0 1, 



2^ 1 y/3 

e 3 = h - — 

2 + 2 ' 
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we find that the characteristic determinant (|1.5[) is given by 



j=0 

2 



■ ip(u> — ui 2 ) 



,r=0 
2 



(3.3) 



in particular A°(p) = ip(ui — ui 2 ) ^^a/e 1 " p . (3-4) 

r=0 

In all these cases, the PDE discrete spectrum is equal to the discrete spectrum of the 
operator (jSmithLlioTlh . 

A calculation of the associated polynomials shows that the differential operator is 
Birkhoff regular if /? 7^ 0. On the other hand , the differential operator S° obtained when 
P = is degenerate irregular by Locker's (|2008l ) classification. 

Although the only difference between the coupled and uncoupled operators is the first 
boundary condition, it is expected from the classification result that the operators have very 
different behaviour. This difference is reflected in the spectral behaviour of the two differential 
operators, as is shown in section 13.11 below. The initial-boundary value problems also have 
very different properties. These are discussed in section [3~2l 

3.1 The spectral theory 

In this section we use operator theoretic results to investigate whether the eigenfunctions of 
form a basis. 

The case P 0. It is shown inlSmithl lj201lh that this differential operator is regular, 
hence by the theory of lLockerl lj2000h we conclude that the eigenfunctions form a complete 
system in T>(S). 

The case j3 = 0. Since this differential operator is degenerate irregul ar, Locker's th eory 
does n ot apply. Indeed, th e proof of the following result can be found in ISmithl (|201ll ) and 
also in [Papanicolaou! (|201lh . 



Theorem 3.1. Let S° be the differential operator corresponding to p = 0. Then the eigen- 
functions of S° do not form a basis in T>(S). 

The proof is based on the following steps: 

• The eigenvalues of S° are the cubes of the nonzero zeros of the exponential polynomial 

e H + uje ' + ui e . (3.5) 



The nonzero zeros of expression (|3.5[) may be expressed as complex numbers Ok, took, 
u) 2 Ok for each k G N, where Re(crfc) = and Im(at) > 0. Then ou is given asymptotically 
by 

- ia k = ( k + I) + O (e-^ nk ) as k -> 00. (3.6) 



V3 V 6 

• For each k 6 N, <j>k is an eigenfunction of S° with eigenvalue o\, where 



r = 

The adjoint operator (S )* has eigenvalues — cr| , k G N, and eigenfunctions 

2 

Mx) = Y. e ~ luircTkX { e ~ lu,T+2,Tk - e ~ lu>r+1 ' Tk ) > fceN - ( 3 - 8 ) 
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• Define ^k(x) = il>k{x) / (il)k,4>k) ■ Then there exists some minimal Y £ N such that 
((0fc)JL y , (*fc)^L y ) is a biorthogonal sequence in AC n [0, 1]. Moreover 

{i> k ,4>k) = (-l^e^+s) +0(1) as fc^oo. (3.9) 

• The eigenfunctions have the same norm and it grows at a greater rate than their inner 
product. 

9 o 3\/3eW fc+1 ^ I e^ k \ 

II^H 2 = ll<fal| 2 = ^ (fc+ i) + °(^H asfc ^~- ( 3 - 10 ) 

• Assume Y = 1 (if Y > 1 the biorthogonal sequence ((</>fc)^L y , (?/>fc)^Ly) is not complete). 
Then the projections Qk = ||0fe|| H^fe || are well defined, and 

Hk\\ 2 



\{>Pk, <t>k)\ 

= -. tt + O — ; asfc^oo. (3.11) 

27r(fc + i) ^ k ) 

Hence the biorthogonal sequence is wild. Now the results of iDaviesl (|2007l . Chapter 3) 
show that (4>k)ksn is not a basis in AC"[0, 1]. 

The case /3 — as a limit 

We now consider the uncoupled case as the limit f3 — t of such calculations for the coupled 
operator. The zeros of Apde 13 are given by 

(k - i) 7r + ilog(-/3) + O ( e _V ?'" r ) even, 



(-fc- |)7r + ilog(-/3) + (^e^— J fcodd, 

and the eigenfunctions of S and S* are given by equations (|3.7p and ()3.8|) respectively using 
the new cr*,. After a suitable scaling, the eigenfunctions of the operator and its adjoint form a 
biorthogonal sequence. 

A direct computation shows that for f3 7^ 0, the fastest-growing terms in || cancel out 
so that, for large k, 



t\\ 2 = 0[e^k^ 



bt^ k )=o{e^k-^ 



(This cancellation does not occur in the case /3 — 0.) This causes the projection operators Q n 
to be uniformly bounded in terms of the parameter j3 for /3 £ [— 1 + s, — e] for every < e < 1. 
However, the bound is not uniform as e — > 0. This lack of uniformity is reflected in the 
transition from a regular to a degenerate irregular problem. 

It is useful to compare the positioning of the eigenvalues, a^. Asymptotic estimates as 
well as numerical evidence suggest that for any particular /3 6 (—1,0) the zeros of Apde are 
distributed approximately at the crosses in Figure[T] the solid rays and line segments represent 
the asymptotic locations of the zeros; the dashed lines are dD, the contours of integration in 
the associated initial-boundary value problem. As /3 — > 0~, hence log(— /?) — > — oo, the solid 
rays move further from the origin, leaving the complex plane entirely in the limit, so that the 
solid line segments emanating from the origin extend to infinity. 



3.2 The PDE theory 

We now show, using the Fokas method, that while the initial-boundary value problems is 
well-posed for any value of /3, the solution admits a series representation only if /J 7^ 0, in 
agreement with the operator theory result of the previous section. 
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The case /3 ^ 

It is already well known (Fo kas and Pellonil . 120051 : ISmithL l201ll ) that in this case we have 
the following result. 

Theorem 3.2. The initial-boundary value problem associated with (5' /3 ,i) is well-posed and 
its solution admits a series representation. 

The case /3 = 

Theorem 3.3. The initial-boundary value problem associated with (S°,i) is well-posed but 
the problem (S°, —i) is ill-conditioned. 

Proof. The proof of the well-posedness claim in this statement can be found in lSmitr] (|201ll ). 
However, for this example we now show that the statement '(^ ± (p)/Apde (p) — > as p — > 00 
from within the sets enclosed by T ± ' does not hold, implying that (S°, —1) is ill-conditioned. 
The reduced global relation matrix in this case is given by 

(c 2 {p) c 2 (p)e-'» Cl (p)e-'P \ 
A(p)= \c 2 ( P ) c 2 (p)e—" ctp)^' , 

hence its determinant Apde (p) = A°(p) given by H3.4|) . and the functions 

2 

Ci(p) = ip(u 2 - u) ^2u r qo(uj r p)e l " p , 

r=0 

2 

* / \ • " I r \ ( r + 1 — i(*j r+1 p r + 2 —iu> r+2 p\ 

Ca(p) = iP2^<lo{u p) lu> e p -u e IJ J , 

r = 

Cs(p) = ipj^q^p) (e-^ +2p - e — , 

r = 

Up) = C«(p) = Ce(p) = 0. 

As a = i, the regions of interest are 

-Ej = Ej\ {neighbourhoods of each crfe}, 
^ = {p G C:^ T ^<arg(p)<^}. 
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We consider the particular ratio 

Cs(p) 



P eE 2 . (3.13) 



For p G E 2 , Re(iuj T p) < if and only if r — 2 so we approximate ratio ()3.13|l by its dominant 
terms as p — > 00 from within E 2 , 

(g (p) - q (ojp))e— + q (a; 2 p)( e — " - e^) + o(l) 
(w 2 - uj)e*P + (1 - ijj 2 )e l "P + o(l) 

We expand the integrals from qo in the numerator and multiply the numerator and denominator 
by e~^ p to obtain 



(3.14) 



Let (7?j)j g N be a strictly increasing sequence of positive real numbers such that pj = 
Rje % ~ £ E 2 , Rj is bounded (uniformly in j and k) away from {2j=(k + ~) : fc £ N} and 
i?j — > 00 as j — > 00. Then pj ~ > 00 from within E 2 . We evaluate ratio (|3.14[) at p = pj, 

l / 1 (2ie^(i-)-^ i 'sin(^) - e -^a-) (l - e'^fj q (x) dx + o (e~^ 

(3.15) 

The denominator of ratio (|3.15[) is bounded away from by the definition of Rj and the 
numerator tends to 00 for any nonzero initial datum. □ 

Remark 3.4. In the proof of Theorem 13.31 we use the example of the ratio g~~7^ being 
unbounded as p — > 00 from within E 2 . It may be shown using the same argument that Apde'Ip) 
is unbounded in the same region and that both these ratios are unbounded for p £ E3 using 
Pj = Rje l ~s~ for appropriate choice of (Rj)j^. However the ratio 

Up) C + (p) 



ApdE (p) ApDB (p) 

is bounded in E\ = E + hence it is possible to deform the contours of integration in the upper 
half-plane. This permits a partial series representation of the solution to the initial-boundary 
value problem. 

Remark 3.5. For all /3 6 (— 1, 1) the final time boundary value problem is ill-posed. The 
asymptotic location of the zeros of Apde , along rays wholly contained within {p 6 C : 
Re(-ip 3 ) < 0} means that for nozero initial data the solution exhibits instantaneous blow-up. 
Nevertheless, for all /3 G (— 1, 0) U (0, 1) the final time problem is well-conditioned. In the case 
/3 — 0, the final-time problem becomes ill-conditioned and 5 becomes degenerate irregular 
under Locker's classification. 

When P = ±1, S is self-adjoint and the initial- and final-boundary value problems are 
both well-posed. For |/3| > 1, the final-boundary value problem remains well-posed but the 
initial-boundary value problem becomes ill-posed. Thus the self-adjoint cases represent the 
transitions between well-posedness of the initial- and final-boundary value problems. Anal- 
ogous to the /? = case, in the limit j3 — 00, the initial-boundary value problem becomes 
ill-conditioned, the solution to the final-boundary value problem may not be represented as a 
series and S becomes degenerate irregular. 
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3.3 Comparison 

The explicit computation of the operator norms in section [3721 requires the evaluation of the 
biorthogonal family of eigenfunctions and the precise asymptotics for the corresponding eigen- 
values. 

On the other hand, the integral representation of the solution of the boundary value prob- 
lem can be constructed algorithmically from the given data, without the need for any precise 
asymptotic information about the eigenvalues, e xcept their asymptotic location (always along 
a ray for odd-order problems; see ISmitbL I2012I . Theorem 6.3). This is sufficient for a direct 
analysis of the terms that blow up and prevent deformation of the contour of integration and 
a residue computation around the eigenvalues, thereby precluding a series representation of 
the solution. 

In the example above, the particular term in the integral representation exhibiting this 
blow-up is the term 

J 2zsm f - - 3 J e 2 1 ' 2 dx ~ — e 2 (3.16) 

where the right hand side is obtained by an integration by parts. Note that in particular we 
can choose Rj = ^| (k + |) . 

Comparing this with expression (|3.10p . 

11/ II* IU II 2 3^ett( fc+ s) (e% k \ 

IWI =M = ^(fe+i) + °[ir) asfc ^°°- 

it is evident that the lack of boundedness of the norms of the operators, responsible for the 
lack of the properties of a basis for the eigenfunctions biorthogonal family, is exactly the same 
lack of boundedness in the integrand of the integral representation for the solution of the PDE, 
yielding a barrier to the contour deformation. Indeed, using the notation of Theorem 12. 121 we 
have shown that, for this example, 

Q fe = O sup -jp f-L and 

\l>|e|>0 Apde (W(7fc + e) 

Q(u>a k + e,ip k ) 

SUP A r— -—T = O (HQfell) • 

l»|e|>0 A PDE (Wfffc + e) 

This is a tighter bound on the blowup of ||Q&|| than that obtained in section [5] No examples 
have been found that violate the tighter bound but an example is presented in section [4] for 
which IIQfeU = O(l) while the spectral ratio grows exponentially with k. 



4 3 order pseudoperiodic examples 

In this section we outline the analysis of another class of boundary value problems, depending 
on a real parameter /3, for the linearized Korteweg-de Vries equation. Namely: 

* = x€[0,l], te[0,T], (4.1) 

q(x,Q) = qo{x), x G [0, 1], 

q{0,t) = q(l,t), 9«(0,t) = -^« B (l,t), q xx (0,t) = q xx {l,t), t G [0, T] /3 G R 

where go G T^{S) is a known function. 

For all P ^ 0, these are pseudoperiodic boundary conditions. In the limit as the constant 
P — > 2, the boundary conditions fall into the special class o f pseudoper iodic conditions for 
which the solution cannot be represented as a discrete series (|Smithl . 120121 . Section 5). As in 
Section [3] the spectral properties of this limiting case are very different from the case P 7^ 2. 

In this section we analyse the behaviour of the associated differential operator in the two 
cases. To avoid technicalities, and to concentrate on the P — 2 limit, we assume in what 
follows that p G (2-e,2]. 
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The associated differential operator 

For the real parameter /3 £ (2 — e, 2], we investigate the differential operator S 13 with pseu- 
doperiodic boundary coefficient matrix 

A -1 \ 
A= 1 /3 , 
\0 1 -1/ 

and the associated initial- and final-boundary value problems II ^ and II'' 3 . 

Remark 4.1. The restriction from j3 £ R \ {-1,0, 1/2} to )3 6 (2 - e, 2] is not of any 
consequence other than notational convenience but the cases /3 = — 1, f3 = and /3 = 1/2 
require special treatment. 

Indeed, f3 = 1/2 is equivalent to the final-boundary value problem II' being well-posed but 
with solution lacking a series representation (as II is ill-conditioned) and, as S 1 ^ is the adjoint 
of S 13 for any /3 7^ 0, the below analysis carries over to this case with a relabeling between S 
and S* . 

If /3 = then the boundary condi tions are no l onger pseudo-periodic. A description of 
well-posedness for this case is given in lSmitbl l|2013al ). 

If /3 = — 1 then the operator is periodic hence, from the classical theory, its eigenf unctions 
form a basis and the problems II and IT' are both well-posed. 



4.1 The spectral theory 

In this section we attempt to use operator theoretic results to investigate whether the eigen- 
functions of S 13 form a basis. 
The case /3 < 2 

It is shown bv lSmitbl l|201ll ) that this differential operator is regular, hence by the theory 
of lLockerl ||200Ch ' we conclude that the eigenfunctions form a complete system in T>(S). 

The case /3 = 2 

Since this differential operator is degenerate irregular, Locker's theory does not apply. 
However, in this example we are unable to apply Davies' method to discer n whether the 
eigenfunctions form a basis. The eigenfunctions form a tame (in the sense of lDavied.l2000l ) 
system, which is a necessary but not sufficient condition for a basis. 

Indeed, following the same outline method as in Section [3. 11 we obtain 

• The eigenvalues of 5* 2 are the cubes of the nonzero zeros of the exponential polynomial 

e~ lp + <T iup + e-^ 2p - 3. (4.2) 

The nonzero zeros of expression (|4.2|) may be expressed as complex numbers at, uj/Jk, 
u) 2 Ok for each k £ N, where Re(crfc) = and Im(at) < 0. Then ou is given asymptotically 
by 

1/Tt. — 

• Let 

2 

cj> k {x) = Y, UJ, ' e ' Ujr ° kX {er iura " - e i " r+2 ' Tk -e™^ 1 "* + l) , k £ N. (4.4) 

Then, for each k £ N, cf>k is an eigenfunction of S 2 with eigenvalue a%.. 

• The adjoint operator (S 2 )* has eigenvalues { — cr^ : k £ N}, corresponding to eigenfunc- 
tions 

2 

and there are at most finitely many eigenfunctions of (S^ 2 ')* that are not in the set 
{V>fc : k £ N}. 



io k = ^(k-^]+0 (e-^/A a s k -> 00. (4.3) 
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Let 



* fc (aO = 



(4.6) 



Then there exists a minimal Y £ N such that ((</>fc);*Ly , (^j s )^ =y ) is a biorthogonal 
sequence in ^4C n [0, 1]. Moreover 



v ^_ e t%( fc -^) + O^^AT 1 ) as fc -» oo. 



(4.7) 



• The eigenfunctions have the same norm and it grows at t/ie same rate as their inner 
product. 

Il^ll 2 = ll^fcll 2 = TTTp^r, e t%( fc -i) +Q( e ^fc- 1 ) as ^ oo. (4.8) 
2(fc- ±)tt 

• Then the projection has norm ||</>fc || ||^fc || , which is bounded uniformly in k. From 
this result, it is impossible to determine whether the eigenfunctons form a basis or not. 

4.2 The PDE theory 

As shown in lSmit hi l|20ll Example 5.2). n is ill-posed if and only if ft = 2. Via Proposition ^. 71 
this yields the result that the analysis of section I^TTI could not — the eigenfunctions do not form 
a basis. 

Prop osition 4.2. Let = 4fc7r/\/3 and let p k = Rue^^ . Then, using the notation o ^Smitti 
\20li) . the ratio 



4 2) (Pk) _ (-l) fc (gT(0)-2g T (l))e^/ 2 



A PDE (Pfe) 

Proof. A quick calculation yields 
^pde (Pfc) = i^SRl 

hence 



SRI 



+ 0(e Rk/2 R^ 3 ), as k -> oo. 



2 2 

3/3 + 3 + (/3 - 2) ^ e iujpk + (1 - 2/3) ^ e - 1 " 3 ^ 

3=0 i=0 



Ap 2 Q B {p k )=izV^Rl 



j=o 



The spectral function 



(4.9) 



(4.10) 



Z 

(p k ) = iV3u 2 Rt^u 2j q T (u j Pk) (e^' Pfc - e~^ + lp " - e — J+ V + l) 



3=0 



is independent of /3. 

By the definition of pk, the functions 



,«fc/2 ifl fc V^/2 



grow exponentially with fc, while 



-i" Pfc _ g-^fc g'Pfc — e 



--Rfc/2 c -i-R fc v / 3/2 giwpj, _ e -R k /2 e iR h V3/2 



decay. Hence 



A pde (Pk) = (-l) fc+1 i6VaRfce K * /2 + 0(e- Rk Rl), as fc ^ oo. 
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Also 



V2 2) (Pk) = -iu 2 V3R 2 k ^-g T (pfc)e lupk - u 2 q T {ujp k )e lf>k + u>q T (uJ 2 >fc)e w Pfc ) 

+ 0(e ii ' t/2 ^) 

= -cAV^ f- jf 1 g T ( a; ) e i ^^ +1+I ^ a: - 1 )) dz 

-u, 2 ^ 9T (x)e^ L[a;+1 - I ^ (3; - 1)1 dx + "£ qT(x)e R « ll -* ] dx) 

+ 0{e Rk/2 R 2 k ) 



2 i\f?>R\ 



2q T (l)e Rk 2 2g T (l)e flfc <?t(0)£ 



fl fc (l + iV3) ' ft&(l-n/3) -Rfc 

+0(e Rk R k 2 ) ) j +0(e Rk/2 R 2 k ) 

= -iVSifc (2g T (l) - g T (0)) e Rk + 0(e Rk ). 

Note that gr(l) = g(l,T) = —q(0,T) = — Qt(0), by the first boundary condition. Hence, 
provided we can be sure that gr(0) 7^ 0, 2qr(l) — ?t(0) 7^ 0. □ 

As < arg(pfc) < 7r/3, and R k was chosen to ensure that Apde (pk) is bounded away from 
0, p k G Dl Hence, bv lSmithl l|20ll Theorem 1.1), n is ill-posed. 

The rate of blowup exhibited in Proposition 14.21 is maximal in the sense that for any 
sequence (p k )ken such that |fffc_i| < \p k \ < |fffc| and for any j € {1,2,3}, 

^P^ = 0(e Rk ^ 2 R k 2 ). 

A PDE (Pk) 

The problem n' is well-conditioned for all ft £ (2 — e, 2]. Indeed, for any sequence (pk)ken 
with p k £ fl r and — > 00, we find the asymptotic behaviour: 

rif(Pk) x 
C(|p fc | ). 



4.3 Comparison 

In order to find the asymptotic behaviour of ||Qfc||, the complex calculation outlined in sec- 
tion 14.11 is necessary. However, the result we obtain is that the projection operators are 
uniformly bounded in norm, from which we cannot discern whether the eigenfunctions form a 
basis. 

The calculation required to prove Proposition 14. 21 is relatively simple and from that result, 
via Proposition 12.71 it follows that the eigenfunctions are not a basis. 



Conclusions 

In this paper, we have gathered and summarised old and new results on a newly analysed 
correspondence between the spectral theory of linear differential operators with constant coef- 
ficients and the analysis and solution of IBVPs for linear constant coefficient evolution PDEs. 
We also presented two specific examples to illustrate the power of this connection for inferring 
results on the spectral structure of the operator. 

In Section O we developed a new method for showing that the eigenfunctions of certain 
differential operators do not form a basis. This method relies crucially upon finding a well- 
posed IBVP whose solution cannot be represented as a series. 

In SectionslSHU we compare the new method to the established method of Davies by apply- 
ing each method to examples. The calculations we present suggest that the PDE approach is 
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more straightforward in deriving estimates for the boundedness of projector operators, hence 
results on the existence of eigenfunction bases. Indeed, it is sufficient to estimate the bound- 
edness of functions constructed algorithmically in certain well defined complex directions. 

The second example represents a case where the new method yields a result but the 
operator-theoretic methods we considered do not. Indeed we show that the solution of the only 
well-posed initial-boundary value problem cannot be represented as a discrete series hence, by 
Proposition 12.71 the eigenfunctions cannot form a basis. But the eigenfunctions are not wild, 
indeed the associated projection operators are uniformly bounded in norm, so we cannot reach 
the same conclusion using e.g. the operator-theoretic framework of Davies. 

The remainder of Section [2] investigates the relation between the two methods. Indeed, for 
the class of operators we discuss, determining the wildness of the eigenfunctions is equivalent 
to the calculation of precisely the same quantities used to determine well-posedness of the 
associated initial-boundary value problems. 

It is expected that the well-posedness of both the initial- and final-boundary value problems 
is sufficient to guarantee that the projection operators are uniformly bounded in norm. 

The applicability of the new method has only been shown for eigenfunctions of the class 
of differential operators considered herein, whereas Davies' method could be applied to any 
complete biorthogonal system, whether it is constructed from the eigenfunctions of differential 
operator or not. How ever, it should be possible to extend the new method, along with the 
results of ISmit 2 <|2012h to a wider class of differential operators, providing a powerful tool to 
investigate the spectral properties of linear differential operators. For example, throughout 
this work we have assumed that S = (—idx) n . A general constant-coefficient linear differential 
operator may have more terms, but its principal part could always be represented by such an 
operator S. As the spectral behaviour of the operator is governed by its principal part, we 
expect the above results to carry over to such operators. 
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A Appendix 

The statements of results I2TB112 . 1 2 1 all require the following additional conditions: 
Condition A.l. The boundary coefficient matrix A is non-Robin: 

None of the boundary conditions represent couplings between different orders of boundary 
function. 

That is, for each k £ {1, 2, . . . ,n}, if auj 7^ or f3kj 7^ then atr = = fikr for all r 7^ j. 

Note the following contrast with Robin's original definition. Our Robin/non- Robin classifi- 
cation is independent of coupling between the two ends of the interval; the boundary condition 
q x {0, t) = 5(1, t) is of Robin type and couples the ends of the interval. 

Condition A. 2. Recall that A is reduced row-echelon form. The boundary conditions are 
such that if the boundary function of order r at one end corresponds to a pivoting entry in the 
boundary coefficient matrix A then the boundary function of order n — 1 — r at the other end 
must correspond to a non-pivoting entry in A. Further, the coupling constants for coupled 
boundary conditions of order r and n — 1 — r are equal. 

For simple boundary conditions, this means that if the boundary function of order r at 
one end is specified then the boundary function of order n — 1 — r at the other end must not 
be specified. 
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